Superconductor- Insulator transition in a network of 2d percolation 
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Abstract - In this paper we characterize the superconductor-insulator phase transition on a 
network of 2d percolation clusters. Sufficiently close to the percolation threshold, for p ~ p c , 
this network has a broad degree distribution, and at p = p c the degree distribution becomes 
scale-free. We study the Transverse Ising Model on this complex topology in order to character- 
ize the superconductor-insulator transition in a network formed by 2d percolation clusters of a 
superconductor material. We show, by a mean- field treatment, that the critical temperature of 
superconductivity depends on the maximal eigenvalue A of the adjacency matrix of the network. 
At the percolation threshold, p = p c , we find that the maximal eigenvalue A of the adjacency 
matrix of the network of 2d percolation clusters has a maximum. In correspondence of this max- 
imum the superconducting critical temperature T c is enhanced. These results suggest the design 
of new superconducting granular materials with enhanced critical temperature. 



Introduction. — Complex topologies strongly affect 
the phase diagram of classical phase transitions [UHJ- 
Scale-free networks, with power-law degree distribution 
P(k) ~ fc~ A and diverging second moment of the degree 
distribution, i.e. A e (2,3], have a phase diagram of the 
Ising model, or of the percolation phase transition, which 
changes significantly with respect to the phase diagram 
of the same models on Poisson networks OH]. Moreover, 
the spectral properties of the networks drive a number of 
other critical phenomena [3H5]. 

Recently, quantum phase transitions defined on complex 
networks are starting to gain a growing attention. Quan- 
tum critical phenomena depend on the topology of the 
underlying lattice, as it has been shown for Bosc-Einstein 
condensation in heterogeneous networks [lOj , for Anderson 
localization on scale-free networks with increasing cluster 
coefficient [HUH], and for the Bosc-Hubbard model on 
complex scale-free networks [13] . Several papers have also 
characterized quantum processes on Apollonian networks 
|14[I15|. which provide an example of scale- free networks 
embedded in two dimensions [16H18] . In this context, the 
Transverse Ising Model is attracting increasing attention. 

The random version of this model, the Random Trans- 
verse Ising model, has been proposed to study the 
superconductor-insulator phase transition in granular ma- 



terials [T3H2"Tj . In each grain of granular materials the 
superconducting order parameter is well denned, and the 
grains are coupled to each other by the pair transfer term 
and the disorder is modulated by different on-site ener- 
gies. Therefore, the physics is similar to the superconduc- 
tivity in Joscphson junction arrays [35]. In [TSU^O] this 
model has been studied on a quenched Cayley tree net- 
work by using the quantum cavity method. Subsequently, 
this model has been studied on annealed complex networks 
|21j . i.e. networks that dynamically rewire their links. It 
has been shown that the phase diagram is strongly af- 
fected by a scale-free network topology of the underlying 
networks on which the model is defined. In particular 
when the second moment of the degree distribution (k 2 ) 
diverges with the network size, the critical temperature for 
the superconductor-insulator of this model diverges. This 
suggests that, by modulating the topology of the underly- 
ing network, the critical superconducting temperature can 
be enhanced. 

Nevertheless, the physical realization of complex net- 
works as the underlying structure where condensed matter 
processes take place, poses many questions. Here, we pro- 
pose to model a 2d percolation pattern of superconducting 
materials as a network formed by interacting supercon- 
ducting clusters. This model is inspired by recent exper- 
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imental data on the structure of high temperature super- 
conductors (HTS). Cuprates, diborides, and iron based 
compounds are made of 2d superconducting layers, inter- 
calated by spacer layers, with a lattice misfit strain [53,24 
that induces incommensurate phases [35] • Here, defects 
self organization [23H7] plays a key role in the formation 
of complex nanoscale electronic heterogeneity [35] . A per- 
colating network of 2d superconducting grains has been 
observed in doped diborides [29], in the electron doped 
iron chalcogcnidcs 1301131] , in cuprates using scanning nano 
x-ray diffraction |32 | l33 j . and time resolved x-ray diffrac- 
tion [31[3S]. 

These data show that high temperature superconduc- 
tivity (HTS) is also controlled by the spatial nanoscale 
arrangement of superconducting grains. Therefore the 
theoretical focus for the mechanism of high T c is shifting 
toward the search for the mechanism enhancing T c in a 
network of superconducting grains in a 2d lattice [3"M3"5] . 

In this scenario, the pseudogap phase, appearing in the 
underdoped phase of cuprates at temperatures T* above 
the superconducting phase, is assigned to disconnected su- 
perconducting clusters formed at T*, which become phase 
coherent at T c [33]. The in plane resistivity of the nor- 
mal state (well below T c in the presence of very intense 
magnetic fields used to destroy the superconducting state) 
exhibits a metal to insulator (M-I) crossover [40] ST] at 
the same doping as the maximum of T c . This transition 
from an insulating to a metallic normal state is assigned 
to nanoscopic superconducting grains in the underdoped 
phase, percolating where the critical temperature for su- 
perconductivity reaches its maximum. Therefore, we ex- 
tend our recent theoretical discovery of the increasing T c 
in scale free networks [2T] to the case of networks of su- 
perconducting grains in a 2d lattice in the proximity of 
the percolation threshold to be compared with experimen- 
tal realizations. In particular, we construct a network of 
2d percolation clusters which are in close proximity. We 
show that at the percolation transition the network be- 
tween clusters becomes scale-free. Moreover, we show, by 
a mean-field calculation, that the critical temperature of 
the Quantum Transverse Ising Model, proposed to study 
the supercondutor-insulator transition in this system, has 
a maximum at the percolation threshold p = p c for any 
given size of the 2d dimensional array. This provides a 
new roadmap for the design of complex superconducting 
materials with enhanced critical temperature T c . 

Network of 2d percolation clusters. In order to 
mimic the fractal background present in cuprates [3"2"ll33j . 
we consider a granular 2d material formed by percolating 
clusters of superconductors that could be realized in artifi- 
cial arrays of superconducting grains connoted by Joseph- 
son Junctions on a 2D surface. We consider a square lat- 
tice of size L where each site is occupied by a grain of 
superconducting material with probability p. The array 
of grains displays a site percolation pattern [42] . Two oc- 
cupied sites, m and n at distance d mn < y/2, belong to 
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Fig. 1: (Color online) Percolating clusters in an array of size 
L — 10 with p — 0.30. The number indicates the labelling 
of each to the five clusters. In a coarse-grained picture any 
two percolating clusters i and j at distance dij < VE are con- 
nected together. Therefore cluster 1 is connected to all the 
other clusters and no other links are present in this network. 
The adjacency matrix of the network is given by Eq. Q. 



the same percolating cluster. In other words a cluster is 
formed by occupied nearest neighbors and next nearest 
neighbors. The percolation threshold for this case was 
studied in [43] , where the percolation threshold was found 
to be p c = 0.407.... 

Here, we define the distance between two clusters, i and 
j, as the minimal distance between any site n belonging 
to the first cluster and any site m belonging to the second 
cluster, i.e. 



ieC(j)dnm 



(1) 



where C(i) are the set of sites belonging to the cluster i. 
In order to have a coarse grained view of this system we 
construct a network of percolating clusters. Each perco- 
lating cluster i is linked to any other percolating cluster 
j which is in close proximity. There are different ways to 
implement this definition. Here we consider two percolat- 
ing clusters, i and j linked, if their distance satisfy the 
following requirements: 



2 < dij < VI. 



(2) 



In Figure [TJ we show a small array of size L — 10 and the 
percolation clusters obtained in the method mentioned for 
p = 0.30. In this figure we have labelled the five perco- 
lation clusters with an integer number from 1 to 5. Each 
pair of percolating clusters i and j are linked in a network 
if the distance between the cluster satisfies the condition 
given in Eq. ([2]). The adjacency matrix a of the coarse- 
grained percolation network, that models the interaction 
of the clusters in Figure [TJ is given by 
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Fig. 2: (Color online) The maximal eigenvalue of the network 
as a function of p for different linear sizes L of the array. The 
data shown are generated from 2d array of linear size L — 
100 and the data are averaged over 1000 realizations of these 
networks. The transition for this array of linear size L = 100 
is about p c = 0.407 .... 



At the percolation transition it is well known that the 
percolating clusters becomes fractal [32]. Moreover it is 
well known that the distribution P(s) of the size s of the 
percolating clusters will scale like 



P(s) ~ s- r $(s/.s c ) 



(4) 



with s c = \p— Vc\~° [42| . This distribution becomes a pure 
power- law distribution at p = p c . As we change the pa- 
rameter p, we can study the degree distribution of coarse 
grained network. What we observe is that, as we approach 
the percolation transition, the degree distribution of the 
coarse grained network becomes scale-free. This is intu- 
itively explained by the fact that larger clusters will have 
a larger perimeter, and therefore more possibilities to link 
with nearby percolating clusters. In Figure[2jwe show the 
degree distribution of the network for p < p c and p > p c 
respectively For p < p c the degree distribution displays 
an exponential cutoff. For p > p c the degree distribution is 
dominated by hub clusters with a very high degree present 
in the network. Finally for p = p c the degree distribution 
is scale-free. In Figure [3] we show the network of 2d per- 
colating clusters as a function of p for a value of L = 100. 
For p < p c the network has many clusters with relatively 
low connectivity, at p ~ p c a hierarchy of hub nodes ap- 
pear, finally for p > p c the number of cluster is strongly 
reduced and one hub node dominates the structure of the 
network. 




Fig. 3: (Color online) The network between 2d percolating 
clusters. From the top to the bottom we show typical networks 
for p = 0.32, 0.40, 0.48. 
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Fig. 4: (Color online) The spectrum of a single network formed 
by a 2d percolation cluster of and array of linear size L = 500. 
In particular we plot the rank iV> (A) of the eigenvalue A against 
A. The rank of the eigenvalues is proportional to the cumulative 
distribution N > (X) oc -P>(A) of the eigenvalues A. We show 
that, at the percolation threshold, when p ~ p c , the rank of 
the eigenvalues scale as a power-law, i.e. N> (A) oc A _a+1 . This 
implies that the probability density of the eigenvalue goes like 
p(A) oc A - a with a ~= 3.1 

Mean field calculation of the Transverse Ising 
model the quenched networks of 2d perco- 
lating clusters. — In this section wc propose to 
study the Transverse Ising model to characterize the 
superconductor-insulator phase transition in granular su- 
perconductors close to the percolation threshold. Wc will 
adopt a coarse grained view of the 2d pcrcolative clusters 
in which each cluster of superconductor material is a node 
of a network, and there is a coupling between nearby per- 
colating clusters. The fact that the degree distribution 
of the network of 2d percolating clusters is scale-free at 
the percolation transition is suggestive, since it was shown 
that on annealed scale-free network the critical tempera- 
ture of the superconductor-insulator phase transition is 
diverging with the network size [3T]. Nevertheless, the 
network under consideration here is both quenched and 
non-random, and we need to investigate further the topol- 
ogy of the model in order to make some conclusions on 
the critical temperature of the superconductor-insulator 
phase transition defined on this network. 

We consider a system of spin variables of , of , for i = 
1, . . . , N, defined on the nodes of a given quenched network 
(the 2d percolating clusters) with adjacency matrix a. The 
Traverse Ising Model is defined as 

" Y.'^- ( 5 ) 

ij i 

This Hamiltonian is a simplification respect to the XY 
model Hamiltonian proposed by Ma and Lee [44] to de- 
scribe the superconducting-insulator phase transition but 
to the leading order the mean-field equations for the order 
parameter are the same, as widely discussed in [T31I2T)] . 
The Hamiltonian describes the superconducting-insulator 
phase transition as a ferromagnetic spin 1/2 spin system 
in a transverse field. We propose to use this Hamiltonian 



to describe, in a granular superconductor, the transition 
from a phase of superconducting grains with no phase co- 
herence (called insulator for granular superconductors) to 
the low temperature superconducting phase with phase 
coherence. Each node of the network corresponds to a 2d 
percolation cluster and the adjacency matrix is the adja- 
cency matrix of the network of 2d percolation clusters de- 
scribed in the previous section. In this model we consider 
clusters with a minimun size large enough so that there 
are no large charging effects and the coherence lenght is 
always smaller than their size. For superconductors with 
short coherence lenght the minimum size could be about 
10-20 nanometers. The spins <Ji in Eq. ([5]) indicate oc- 
cupied or unoccupied states by Cooper pairs or localized 
pairs; the parameter J indicates the couplings between 
neighboring spins, e is the on-site energy. Finally, in this 
model the superconducting phase corresponds to the ex- 
istence of a spontaneous magnetization in the x direction. 
We could also consider a different version of the model 
in which the network of interaction is weighted, and the 
link's weight between two nodes depends on the number of 
sites in the two clusters that are at distance € [2, v5]. 
We have checked that this modified version of the model 
doesn't change the main result of the paper. Therefore, 
for simplicity, in this paper we consider only simple net- 
works with adjacency matrix elements a^- = 0, 1. Close to 
the percolation threshold, the spectrum of the adjacency 
matrix of the network of 2d percolation clusters develops 
a power-law tail p(X) oc X~ a (See figure [4}, therefore los- 
ing memory of the underlying 2d dimensional structure. 
Therefore in order to study the dynamics of the granular 
superconductor, wc perform a mean-field approximation 
in which we put 

of o* c of (oj) + (of >o* - (of) (o*). (6) 

We can therefore consider the following mean-field Hamil- 
tonian, in which we use mf = (of) 

h mf = -jJ2 "■J rr - 1 "' "<./'"-'< - e& i- ( 7 ) 

ij ij i 

In the mean field approach the partition function for this 
problem is given by 

Z = Tr e - pkMF (8) 

with the Hamiltonian given by Eq. (J7J. The mean- field 
Hamiltonian can be written as 

H M f = E + - ^ "■!'"' >"j (9) 

ij 

with 

E = -jY t ofh? F -hY t o?, (10) 

i i 

where the local fields hf IF are given by 

3 
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Fig. 5: (Color online) The maximal eigenvalue A of the network 
as a function of p for different linear sizes L of the array. Data 
are averaged over 500 realizations. 



For this problem the magnetizations along the axis x, mf , 
and along the axis z, mf , can be calculated by evaluating 



Trafe-^' 
Z 



Performing these calculations we get 
Jhf F 



(12) 



tanh (l3^{Jhf F ) 2 +eA , (13) 



^J{Jhf 1F y + e 2 



tanh [f3J(Jhf F ) 2 +e 2 ) . (14) 



The Eqs. ([131 are the self-consistent equations that to- 
gether with Eqs. (| 1 lj) determine the local magnetizations 
mf. Close to the phase transition mf <C 1 the self- 
consistent equation for the magnetization mf is given by 



if = -/^T </,,»': 



x tanh(/3e) 



(15) 



If we diagonalize this equation along the eigenvalues of 
the adjacency matrix a, we find that the phase transition 
occurs at 



1 = J A 



tanh(/3e) 



(16) 



where A is the maximal eigenvalue of the adjacency matrix 
a of the network. 

Therefore for A->oo with N — >• oo then /3 — > for any 
fixed value of the coupling J > 0, and the critical temper- 
ature for the paramagnetic ferromagnetic phase transition 
T c diverges. 

Superconductor-insulator phase transition on a 
network of 2d percolation clusters. In order to 
study the superconducting insulator phase transition on 
granular materials close to percolation, we consider the 
network formed by 2d percolating clusters. The criti- 
cal temperature T c of the superconductor-insulator phase 




Fig. 6: (Color online) Critical temperature T c of the 
superconductor-insulator phase transition given by Eq. (|16|l 
with J — 1 and e = 1 plotted as a function of p. The data are 
averaged over 500 realizations. 



transition provided by the mean-field calculation Eq. (|T6[) . 
depends on the maximal eigenvalue of the adjacency ma- 
trix A of the network. Therefore, we have studied nu- 
merically the maximal eigenvalue of the adjacency ma- 
trix of the network of 2d percolating clusters. In Figure 
[5] we show the maximal eigenvalue A as a function of p, 
and we observe that this eigenvalue has a maximum at 
p = p c = 0.407. . . for L = 100. Moreover, in this figure 
we show that in the network of 2d percolating clusters, 
sufficiently close to p c , the maximal eigenvalue of the ma- 
trix A — > oo as N — ¥ oo. Moreover, for every fixed value 
of N, we have that the maximal eigenvalue A has a peak 
for p = p c . In Figure |6] we show the value of the critical 
temperature given by Eq. (|16|) for the superconductor- 
insulator phase transition on the network of 2d percola- 
tion clusters as a function of p showing a peak for p = p c . 
Therefore we predict that the system will display a max- 
imum of the superconductor critical temperature at per- 
colation. 

Conclusions. — In this paper we have investigated 
how a granular material close to percolation can be de- 
scribed at the coarse-grained level as a complex network. 
We have shown that this complex network can be con- 
structed by linking together the next-nearest clusters. 
Close to the percolation threshold, the distribution of the 
sizes of the cluster is power-law distributed and also the 
degree distribution of the coarse-grained network is scale- 
free. Moreover, the adjacency matrix of the 2d percola- 
tion network has a maximal eigenvalue that diverges with 
the network size close to the percolation transition. In 
this coarse-grained network picture of interaction between 
nearby percolation clusters, we have studied the mean- 
field of a Transverse Ising model. This model is proposed 
in order to characterize the superconductor-insulator tran- 
sition in a 2d array close to percolation. We have shown 
that if the maximal eigenvalue of the adjacency matrix of 
the network diverges, the critical temperature of the sys- 
tem diverges as well. Therefore with this coarse-grained 
network model we predict that granular array show a max- 
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imum of the critical temperature for the superconductor- 
insulator phase transition at the percolation threshold. 

The present theory can now be tested producing artifi- 
cial geometrical arrays of 2d superconducting grains on a 
flat surface. In fact, artificial arrays made of mesoscopic 
superconducting islands placed on normal metal surface 
can now be made by using available advanced nanotech- 
nology [45]. These artificial arrays today provide practi- 
cal tunable realizations of two-dimensional granular (2d) 
superconductivity. An alternative approach is the pro- 
duction of artificial superconducting puddles embedded 
in a 2d layer intercalated by different block layers in a 
layered material. The manipulation of the arrays of su- 
perconducting grains can be done by controlling defects 
self-organization by external fields [34l[35]. Finally these 
results open new perspectives for materials scientists look- 
ing for new room temperature superconductors by focus- 
ing on materials design by pattering at the nanoscale. 
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